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ABSTRACT. We present a version of the domino shuffling algorithm (due to Elkies, Ku- 
perberg, Larsen and Propp) which works on a different lattice: the hexagonal lattice su- 
perimposed on its dual graph. We use our algorithm to count perfect matchings on a 
family of finite subgraphs of this lattice whose boundary conditions are compatible with 
our algorithm. In particular, we re-prove an enumerative theorem of Ciucu, as well as 
finding a related family of subgraphs which have 2( n+1 ) perfect matchings. We also 
give three-variable generating functions for perfect matchings on both families of graphs, 
which encode certain statistics on the height functions of these graphs. 



1. INTRODUCTION 



The theory of perfect matchings on planar bipartite graphs is quite rich and mature, and 
has seen a great deal of activity over the past two decades. The centrai questions in this 
theory are counting questions, in which one attempts to give a generating function, or even 
just a count, of perfect matchings on a fìxed planar bipartite graph G. 

In the 1960s, Kasteleyn (113] developed a powerful algebraic tool for answering (in prin- 
ciple) ali such questions. In particular, one can compute the number of perfect matchings 
of a planar bipartite graph by taking the determinant of a signed version of its bipartite 
adjacency matrix (usually called the Kasteleyn matrìx). 

This paper is concerned with a different phenomenon which is not immediately ex- 
plained by Kasteleyn's formalism. Given a periodic planar bipartite graph, it is often 
possible to fìnd a family of finite subgraphs for which the number of perfect matchings 
decomposes into small, simple factors (as do the generating functions of the graphs). The 
canonical example of this phenomenon is the Aztec Diamond of order n [8] , which is a 
diamond-shaped region of the square lattice. The Aztec diamond of order n has n(n + 1) /2 



Figure 1 . The dP 3 Lattice 
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Figure 2. Diamonds of orders 3 and 3.5 




perfect matchings [8 1. Elkies et al. prove this by using a technique called domino shuffling, 
which is a type of random mapping on the set of perfect matchings on a square lattice. 

This paper follows the approach of [8|, but on a different lattice: namely, the hexagon 
lattice superimposed upon its dual triangular lattice (see Figure[T]i. We name this lattice the 
dP% lattice, since it is associated to the del Pezzo 3 surface (the nature of this association is 
explained in Section[7]i. This association was first made clear in the physics literature (9l in 
a quite general setting; it is from this paper that we take the terminology dP%. Indeed, (9) 
introduces several quivers which correspond to the dP$ surface; ours is their Model I. 

We describe a type of domino shuffling which works on the dP$ lattice. We also define 
an analogue of the Aztec diamond of order ne |Z for this lattice (see Figure [5]), which is 
well behaved under domino shuffling, and prove that the number of perfect matchings on 
these graphs is a power of 2: 

Theorem 1. Let m e |Z and let n = [m\. Then the number of perfect matchings on a 
diamond of order m e ^Zis 

f 2 n (™ +1 ) ifm e N 

\ 2(™ +1 ) 2 ifm+^e N. 

Moreover, we can give a generating function for these matchings; we defer the statement 
of that theorem until Section[8] due to the number of defìnitions required to state the result. 

We constructed our domino shuffle by reverse engineering a particular application of 
urban renewal, and replacing the averaging procedure (replacing it instead with a careful 
accounting of the locai structure of the perfect matching). 

Propp later refìned his domino shuffle in [21 1 by reworking domino shuffling so that it 
can be applied locally, to a square face of a weighted planar bipartite graph. The resulting 
technique, urban renewal, involves altering the graph at a square face S, as well as aver- 
aging the weights of the edges around S in a certain way. Ciucu gave a further refìnement 
to this formalism, reproving the results in as well as many others. Unbeknownst to us 
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when we began this work, Ciucu had already proved half of our Theorem [T] in [5] using 
a powerful reduction theorem. Moreover, there is a type of domino shuffle that one can 
do on any finite planar bipartite graph G: simply embed G inside a suitably large Aztec 
diamond and use the generalized domino shuffle [21 1 on the square lattice. It is therefore 
important to note that our main contribution is the construction of the domino shuffle on 
the dPa lattice itself, rather than its applications to counting Aztec diamonds. Our shuffle 
is indeed different from the ones used in the aforementioned papers. 

There are at least two good reasons to look for simple domino shuffles on particular 
planar bipartite graphs rather than appealing to general theory. The first of these carne to 
light in E2l l24ll . and involves a curious interplay between the combinatorics of perfect 
matchings and algebraic geometry. By studying a particular class of perfect matchings 
on (infinite) periodic planar bipartite graphs, one can compute a Donaldson-Thomas (DT) 
partition function of the toric surface associated to the graph by the construction in J5]- 
The DT partition function is a virtual count of curves in the space (really a type of Euler 
characteristic on the moduli space of subschemes of the surface). There is an emerging 
picture in which domino shuffling corresponds to a geometrie procedure called wall cross- 
ing 10] [16] [18] |23], and which, for example, links the DT theory of a singular surface to 
that of its crepant resolution. The generating functions that we compute in Section[8]are 
a prelude to our next project, wherein we will compute the Donaldson-Thomas partition 
function for this lattice; this turns out to be a generating function for a different type of 
dimer confìguration on this lattice, but which uses essentially the same weighting. We do 
not address Donaldson-Thomas theory at ali in this article, but see [22, 24] for an account 
of the corresponding problem on the square lattice. 

The second modem setting for domino shuffling is in the work of |fl9l [2] and others, 
which describe square-lattice domino shuffling as a type of random point process. This 
process turns out to be partially determinantal, which makes it one of a select few three- 
dimensional statistical physical models which can be treated completely mathematically, 
without applying to heuristic physical reasoning. One current frontier of exactly solvable 
statistical mechanics is precisely this jump from two to three dimensions, and it is quite 
important to understand fully the phenomenology of any tractable models in this area. It 
seems very likely that our domino shuffle fìts into the framework of |2|, but we have yet to 
exploit this connection. 

In a very recent string theory paper iflOl . Franco, Eager and Romo have reproduced 
(up to a change of variables) and substantially generalized our Theorem [T] using different 
methods. We gratefully acknowledge their help in spotting a typo in Theorem [T] in an 
earlier version of this preprint. 

2. The Lattice 

Before taking the time to describe these subgraphs, we must better understand the dP$ 
lattice itself (see Figure [TJ. It is periodic, bipartite, and planar, where each face is bounded 
by a cycle of size 4. The edges incident to centers of hexagons have length 1, and other 
edges have length \/3/3. Color the vertices of degree four black, and the others white. 

Deflnition 1. We define a square to be a face bounded by a cycle of size 4. A long edge 
of a square is one of length 1; a short edge has length V3/3. Every square has a unique 
vertex of degree 3. The tail afa given square is the edge which is incident to this unique 
vertex and not on the boundary of the square. It is incident to another vertex which is not 
the center of a hexagon. See Figure^ 
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Figure 3. Short edge, long edge, tail, square and kite; orientations. 




Figure 4. A strip 




We orient our lattice such that for each hexagon, a white vertex is directly north of the 
center of the hexagon. The square which is bounded by both the center and this white 
vertex has orientation N, and the square directly south has orientation S. We shall often 
have cause to pair up orientations which are in opposite directions on the compass: N/S, 
NE/SW, NW/SE, and we shall speak of these as pairs of orientations . 

Definition 2. A kite is a square together its tail, and a strip of kites is a set of adjacent 
kìtes with the sartie pair of orientations, which meet at a degree-four vertex oftheir squares 
(not their tail). Given a kite K, let the unique vertex ofdegree 3 that is on the boundary of 
the square and incident to the tail be K's root vertex. See Figure^ 

Note that if a kite has its square oriented N/S, and the tail is north of the square, we have 
a N kite, while we have a S kite if the tail is south of the square; the sanie holds for NE, 
SW,NW, SE kites. 

Definition 3. A matching is a set of vertex-disjoint edges of a graph G. It is perfect when 
it covers every vertex of G. 

Two graph transformations which are often used in perfect matching enumeration are 
urban renewal and double-edge contraction (see Figure [5}. Double-edge contraction is 
the simpler of the two techniques, in which a vertex of degree two is merged with its 
two neighbors; it is easy to check that this transformation induces a bijection on perfect 
matchings. 

Urban renewal is slightly more complicated; it replaces a square S in the graph with a 
smaller square 5", together with four new edges connecting the vertices of S to the vertices 
of S' (see Figure |5j. When applied to an edge-weighted graph, it is possible to defìne new 
weights in such a way that urban renewal leaves the perfect matching generating function 
invariant up to a Constant [21]. The reader should note that in our paper, we do not do this 
change of edge weights, and indeed we are often not working with weighted graphs at ali; 
for us, "urban renewal" simply means the graph transformation shown in Figure[5]l. 



DOMINO SHUFFLING FOR THE DEL PEZZO 3 LATTICE 

Figure 5 . Urban renewal and doublé edge contraction 
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Figure 6. The diamond of order 1, D x 




3. The Integer-order Diamonds 

We now define a family of subgraphs called diamonds. This family is indexed by 
theorder of the diamonds, which will be an element of ^Z>o. Our shuffiing algorithm 
will shuffie the edges on a perfect matching on an order n diamond, and outputs a perfect 
matching on an order- (n + 1) diamond. 

We begin by discussing the integer-order diamonds, which were first discussed by 
Propp [20 1 and further studied by Ciucu [5 1, under the name Aztec dragons. 

Definition 4. A diamond of order 1 is denoted Di (see Figure^; it is a connected sub- 
graph of dP^ which consists of 3 squares. One of them is N, while the other two are 
oriented NW and SE and are adjacent on the NE/SW strip they He on. Each diamond con- 
tains a unique degree-four vertex which is on the boundary of ali of its squares. We cali 
this vertex the center vertex of a D\. 

Note that a N-S strip is tiled by copies of Di, half of which have been rotated 180 
degrees. We refer to these as North and South Dis, according to the direction of their 
N-S kite. 

Now, we develop a coordinate system such that we can define a diamond of order n, 
n E N. Recali that the dP% lattice decomposes into strips, each of which can be tiled by 
copies of Di. 

Definition 5. First, fix in the lattice a Di and cali it T(0, 0). If n is even, let T(0, 0) be 
oriented South, and ifn is odd, let T(0, 0) be oriented North. Cali T(l, 0) the Di directly 
east o/T(0, 0), and T{— 1, 0) the Di directly west o/T(0, 0). Similarly, the Di directly 
north o/T(0, 0) is T(l, 0) and the one directly South o/T(0, 0) is T(— 1, 0). If i, j G Z, 
then T(i,j) must be directly north ofT(i, j —1), directly south ofT(i, j + 1), directly east 
ofT(i — 1, j), and directly west ofT(i + 1, j). 
For i € N, let Si be the strip containing T(0, i). 
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Finally, let 



D n = 



U 



-n<i,j <n 



l»l+lil<r 

where the union denotes the union ofvertex and edge sets of subgraphs ofdP^. 



Figure 7. Coordinates on an integer- and half-integer-order diamonds. 
In both, T(2, 1) is shaded. Not ali graph edges are shown. 
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4. The Half-integer-order Diamonds 

We construct a second, new family of finite subgraphs of dP^. These are half-integer- 
order diamonds, D n+ i where neH. 

Again we describe them with a coordinate system. 

Deflnition 6. Let T(0, 0) be oriented North when n is even, and let it be oriented South 
when n is odd. Let D[ denote a Di which has been reflected in the y axis. Decompose 
the N-S strips into copies of D[; let T(i,j) denote the which is i steps east and j steps 
north o/T(0, 0) as in Definition^ 

Define L n to be the SE-oriented square ofT(n, 0), and R n to be the SW -oriented square 
ofT(n + 1, 0). Let D n+ i be the following graph (see figure^: 



( 



D, 



\ ( 



U T dJ) 



— n<i,j<n 

\ l*l+b'l<™ 



U T &3) 



I 0<i<n 
\0<i<n+l+j 



U R n U L n . 



ì 



Note that these graphs no longer decompose completely into copies of D[: there are 
two extra squares, L n and R n . Finally, we specifically describe Di. 



Deflnition 7. Di consists of a single NE square. 
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Figure 8. Applying urban renewal and double-edge contraction to a 
kite 



a 




b 
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Definition 8. Let M n be the number of edges in a perfect matching on a diamond of arder 



Lemma 2. Let n = [m\ . Then 

_( n(3n + l) meN 
Mrn ~\ 3n 2 + An + 2 (m + |) G N. 

Proof. The number of edges in a perfect matching is half of the number of vertices in the 
graph, so this is a straightforward computation using Euler's formula for planar graphs and 
the definì tions of the graphs D n . □ 



Throughout this section, we let n = \_m\ . 

We define a procedure called domino shuffling, whose input is a perfect matching M on 
a D m , and whose output output is a perfect matching on D m , i . We were inspired to work 
out this algorithm by the following observation: 

Lemma 3. Urban renewal applied to ali kites of a given pair of orientatìons, followed by 
contracting doublé edges wherever possible, simply translates the dP% lattice. 

Proof. Without loss of generality, let us perform urban renewal on North and South kites; 
the same argument applies to NE/SW and NW/SE kites. 

First of ali, we choose a N/S kite and cali it K; its square has vertices a, b, c, d and its 
tail has edge de. Then, apply urban renewal to ali of the N/S kites (see Figure[8]l. For each 
kite we have a new square a' ,b' ,c' , d! and 4 new edges aa! , bb' , ce' , dd' . 

We consider the new edges associated with K. We observe that one of the new edges, 
dd' , is adjacent to the root vertex of K . Two other new edges, aa' and ce', are adjacent 
to the new edges of the N/S kites which lie on the same N/S strip as K. We now apply 
double-edge contraction to the vertices a, e, d. We do this to every kite we applied urban 
renewal to. Hence, for each kite we have deleted 1 old and 3 new edges, and added 1 new 
edge. This new edge happens to be bb', which is adjacent to a vertex opposite the originai 
root vertex. Therefore, the new figure is in fact a kite with tail bb', which is oriented in 
the opposite direction of the originai kite. We have fiipped the kite. Moreover, this flip has 
occurred on every N — S kite, hence we have shifted the lattice. □ 



n, such that n e N or (n + |) e N. 



5. The Rules of the Shuffling Algorithm 



8 CYNDIE COTTRELL AND BENJAMIN YOUNG 



Figure 9. Applying urban renewal and double-edge contraction to N\S 
kites in a strip. 




Figure 10. Steps 1 through 3 of shuffling 




Now, we proceed to define the algori thm as applied to a perfect matching M on D m . 
Let n = \m\ . 

Algorithm 1 (The domino shuffle). Begin with a perfect matching on D m . 

(1) If m is an integer, draw an ovai around every NE — SW kite in D m , including 
partial kites. Otherwise, draw an ovai around every NW — SE kite, including 
partial kites. 

(2) Fili in ali of the partial kites around the boundary, by adding new edges to the 
diamond. None of these new edges should be added to the matching, with the 
following exception: the mot vertex of each circled kite is matched. Ifnot, add the 
tail ofthat kite to M. We will add In + 1 tails, by lemma 6 below. 

(3) Flip each kite in each ovai so that it has an apposite orientation. 

(4) Draw the new matching M' using the rules in Algorithm^( see below) and shown 
in Figure \TT\ 

(5) Drop the unmatched vertices on the boundary ovals ( the tail on every other kite is 
never part of the matching.) 

Observe that the edges not on K are not altered by our algorithm: they remain as they 
are (either matched or unmatched). 

As stated earlier, the shuffling is inspired by applying urban renewal and doublé con- 
traction to the squares of the circled kites (but not changing any of the weights on these 
edges - indeed, ali of this will work on an unweighted graph), and keeping track of what 
happens to the associated perfect matching. The rules are shown graphically in Figure 1 1 
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FIGURE 1 1 . Steps 4 and 5 of the domino shuffle 

which is perhaps the most useful format. However, to ensure clarity, we will also describe 
the rules more formally. 

We cali our kite K and adopt the notation of Lemma [3] illustrated in Figure [8] the 
vertices of K are a, b, c, d, e and after urban renewal, the square of K is replaced with a 
new graph H with vertices a, a' , b, b', e, d , d, d' . 

Algorithm 2 (Rules for edges). Let K be a kite with edges as described above. 
Case 1 : No edges in K are in M. 

Then, our root vertex is not matched as it is not adjacent to any edges outside, 
so this case does not occur. 
Case 2: 1 edge in K is in M. 

a. Suppose that the tail ed is in M. Take eitìier a'b' and c'd', or or a! d' and db' 
to be in the new matchìng. This is called a creation. 

b. Suppose that the tail is not in M. Then M contains a short edge, eitìier ad or 
de. Take b'd (respectively, a'b') to be in the matching. We cali this the short 
edge case. 

Case 3: 2 edges in K are in M. 

a. Suppose that the tail of K is in M. Then, either bc or ab is in M. Take a'd' 
( respectively, dd') to be in the new matching, as well as the new tail bb'. We 
cali this case the long edge case. 

b. Suppose that the tail of K is not in M. Then take only the new tail bb' to be 
in the new matching. We cali this case an annihilation. 

Let us analyze this lemma, to see that it works as advertised. 

Lemma 4. Shuffling a perfect matching on ali of dP%, choosing one of the two outeomes 
ofeach creation arbitrarily produces another perfect matching on the translation ofdP 3 . 

Proof. We examine one ovai, and proceed by cases as in Algorithm [2] Our strategy is to 
perform the operations of urban renewal and doublé edge contraction, each of which comes 
with an induced map on perfect matchings. In each case, we verify that we have duplicated 
the shuffling rules of Algorithm[2] 
Case 1 : No edges in K are in M. 

As remarked above, this is not possible. 
Case 2: 1 edge in K is in M. 

a. Suppose that the tail ed is in M. Then, a, b, and c are matched outside of K. 
When we apply urban renewal, we must match the vertices on S'. Hence, we 
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can take either a'b' and c'd! or a'd' and db'. There are two choices. When we 
apply double-edge contraction we see that K' will have either a'b' and c'd! 
or a'd' and db' in its matching. 
b. Suppose that the tail is not in M. In order to ensure that our root vertex 
is matched, we only have a case where our matching contains a short edge, 
either ad or de. Without loss of generality, suppose that this edge is ad. Then, 
b, c, and e are matched outside of K. We only need to match S' and a and 
d on H. Hence, we let aa' and dà! be in the matching as well as b'd . Then, 
we have a perfect matching. Finally, we apply double-edge contraction and 
obtain K' such that only b'd is in the matching. 
Case 3: 2 edges in K are in M. 

a. Suppose that the tail of K is in M. Then, we must have a long edge in M. 
Without loss of generality, let this edge be bc. Then, we have that a is matched 
outside of K. We must now match S' and b and c in H. Thus, we take the 
edges a'd' , bb', and ed . Applying double-edge contraction, we get that the 
edges in K' which are in our new matching are a'd' and b'b. 

b. Suppose that the tail of K is not in M. Then, the 2 edges in M must be a 
short edge and a long edge on S. Without loss of generality, we let these two 
edges be ad and bc; the other possibility is choosing de and ab. Now, we have 
that e is matched outside of K, so we include the edges aa', bb', ed, and dd'. 
Then, when we apply double-edge contraction, we will only have bb' in our 
new matching. We see that this occurs no matter which pair of edges on S we 
select. 

□ 

Lemma 5. When applying step 3 of our shuffling algorithm, we add 2n + l tails to ensure 
that every root vertex is matched. 

Proof. Consider first the case n G Z. For each strip, we see that the two Z?is on the 
boundary are either both oriented N, or both oriented S. If they are both N diamonds, we 
circle two kites on the eastern boundary, and one on the western boundary. If they are both 
S diamonds, we circle two kites on the western boundary, and one on the eastern boundary. 
Two of these kites have their tails in D m , while the third one only has one long edge in 
D m . Hence, its root vertex will never be matched. So, we will have to add as many tails as 
the number of strips. 

Finally, note that there are two more kites on the boundary. One is north of T(0, ri — 1)) 
the other is south of T(0, — n + 1). Each of these kites only have one long edge in D m , 
hence their root vertex also must be matched. 

Thus, we add exactly 2n — 1 + 2 = 2n + 1 vertices. 

The same argument holds for half-integer-order diamonds, if we replace Di by D[ 
throughout. □ 

Lemma 6. Given a perfect matching M on D rn , applying shuffling to M gives us a match- 
ing M' which covers ali ofthe vertices on a diamond of order m + |. 

Proof. First of ali, by Lemma|3] domino shuffling applied to the entire lattice behaves as a 
simple shift. Hence we simply need to show that the number of strips on the new subgraph 
is correct, that the correct types of squares are present on the boundary, and that step (5) 
is deterministic and removes only tail vertices from the boundary. We will do in detail the 
case that rn = n are integers; the half-integer case is similar. We encourage the reader to 
refer to fìgures[2]and[5]to verify these claims. 
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First, note that an N square transforms into a SE square and vice versa; a S square 
becomes a NW square and vice versa; and a SW square becomes a NE square and vice 
versa. 

When we flip the kite directly north of T(0, n — 1) we have added a square that is in a 
strip which was previously not included in our subgraph. This new square is a S W, since the 
originai is NE. Similarly, when we flip the kite which is directly south of T(0, — (n — 1)), 
we include a square which was not on any strips that were on our originai graph. The new 
square will be NE, since the originai was SW. We see that we have increased the number 
of strips by 1 . 

Now, we check that our new boundary is correct. We start with the northernmost kite 
that is added in Step 2: it is a NE kite. Hence, the square we obtain after shuffling is a SW 
square. The northernmost square on the originai subgraph was a SE square, and this one 
will become a N square. There is also an adj acent SW square in the originai graph which 
becomes a NE square. These three new squares share a common vertex which becomes 
the center vertex of an eastern N D\\ these form the correct shape for the top corner of a 
half-integer-order diamond. Moreover, one can apply this to the D\ 's on the strips which 
He on or to the north of the center strip Sq. 

We can apply a similar argument for the strips which lie south of S . We start with the 
southernmost square which is added in Step 2: it is a SW kite. Hence, the square we obtain 
after shuffling is a NE square. The southernmost square on the originai subgraph was a 
NW square, and this one will become a S square. Similarly, we added a NE kite in Step 
2; such that its square will be a SW kite after we apply shuffling. These three new squares 
are such that they share a common vertex which becomes the center vertex of a western S 
Di. Hence, we have the exact shape which we were looking for. One can apply this to the 
Dis on the strips which lie south of So. 

Next, we check that step (5), in which unmatched boundary vertices are removed, is 
compie tely deterministic. Indeed, ali tail vertices on the boundary are always unmatched 
after shuffling. This follows from an inspection of the shuffling rules. Let K be a kite. 
Observe that tail is unmatched within K after shuffling if the vertex v incident to the two 
long edges of K is matched within v before shuffling, which is always the case when v is 
on the boundary of D n . □ 



6. Proof of Theorem 1 

Our proof is essentially an analysis of the domino shuffling algorithm. In short, the 
algorithm behaves much like a l-to-2' n ' +1 correspondence for enumerative purposes. 

Proof. Let n — \m\. When we apply the shuffling algorithm to M E PM(D m ), we add 
2n + 1 tails, but we know that the associated matching on D m+ i has 3n + 2 more edges 
than those in M. We can account for 2n + 1 of these new edges by the 2n + 1 tails added 
when applying our algorithm. Hence, the difference between the number of creations and 
number of annihilations is n + 1. 

Suppose that M has k annihilations. There is a class of exactly 2 k tilings which have 
k annihilations on these kites, and are otherwise identical to M . Moreover, M must 
experience k + n+l creations after shuffling. On each kite which had a creation, there are 
2 possible outcomes. Hence, the 2 fe matchings with k annihilations map to exactly 2 k+n+1 
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Figure 12. Case m = \ 
Figure 13. A fundamental domain of the dual quiver of the dP 3 lattice 



matchings on PM(D m+ i). As such, 

\PM{D n+l )\=T +1 \PM{D n+ì )\- 

\PM(D n+k )\=T +1 \PM{D n )\. 

The count of perfect matchings now follows by a standard induction argument, the base 
cases being PM(Dq) = 1, because the empty graph has only the empty matching, and 
PM(Di) = 2, because the order | diamond is a SW square and has two matchings (see 
Figure [Ù]). " □ 



7. Height function 

Our next aim is to give generating functions for the perfect matchings on D n . To de- 
scribe the statistics encoded in our generating functions, we will need to introduce the 
notion of height function. This is a rather general concept associated to any perfect match- 
ing on a planar bipartite graph (see, for example, [14 1), but it is necessary to work out the 
precise nature of the height function on dP^ to proceed further. 

If a planar graph G is bipartite, its dual graph Q inherits a naturai orientation: orient 
each edge of Q so that a black vertex of G is on the left (say). We use the letter Q because 
the dual graph should be regarded as a quiver ll22l [24], but for combinatorists a quiver 
is precisely the same structure as a directed graph. If we do this construction for the dP% 
graph, we obtain an orientation on the (6,4,3) lattice with arrows pointing counterclockwise 



on hexagons and triangles, and clockwise on squares (see Figure 13 1. 
We will assign an integer to each arrow e of G" as follows: 



(1) H(e) 



1 if e crosses a long edge 

2 if e crosses a short edge 
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Figure 14. Tiles for encoding the height function on the dP^ lattice. 
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This integer is called the height change of e. Moreover, given a perfect matching M on G, 
we will define a function h : V(G') — > Z as follows: fix a to be one of the vertices in G" 
and set H(a) = 0. Now, if e : b — » c is an edge of g' and h(b) has been defined, set 



(2) h(c) = 



h(a) + H(e) - 6 if e € M 
h(a) + H(a) if e & M. 

It is easy to check that h is well-defined up to the initial choice of height h(a). This is 
essentially the same height function defined in lfT4Ì . 

Definì tion 9. Let f be a face of G, bounded by edges E = {ei, • • • ,e2k}- Let M be 
a perfect matching on G such that M fi E = {ei, . . . , e2k-i}- The plaquette flip of 
M around / is the perfect matching M' which agrees with M except around f, where 
M'n/ = {e 2 ,...,e 2fc }. 

It is easy to check that performing a plaquette flip at a kite S changes the height at S by 
±6, leaving the height function otherwise unchanged. Observe also that one can transform 
any perfect matching to any other by a sequence of plauqette flips. This is because if hi , hi 
are height functions for two perfect matchings, then so is min(hi, h^)', one construets its 
associated perfect matching by superimposing hi and hi, producing a collection of disjoint 
doubled edges and loops of even length. 

There are two ways to remove every second edge from a loop, leaving a perfect match- 
ing on the vertices of the loop; for each loop, remove the set of edges which minimizes the 
resulting height function inside the loop. Doing this, we obtain the matching associated to 
min(hi, h-ì). 

In particular, the function 

min h(M) 

MePM n 

is a height function. Observe that this matching admits only plaquette moves which in- 
crease h. 



14 



CYNDIE COTTRELL AND BENJAMIN YOUNG 



In light of these observations, we may visualize a perfect matching on D n as a pile 
of special bricks, contained in a V-shaped tray (just as is done in |8|). Plaquette flipping 
corresponds to adding a brick, and the tray encodes the shape of the minimal height func- 



tion(see Figure 15 See Figure 14 for a picture of the tiles. They are constructed so as to 
meet the following criteria: 

(1) Each tile covers exactly the four faces in Q which are dual to the vertices compos- 
ing a square in G; 

(2) If one tile is placed partially atop another, covering two adjacent kites a, b of G 
joined by a dual arrow e in Q, then the vellicai distance between the two tiles is 
the height change h(e); 

(3) The tiles, viewed from above, have two dimers stencilled on top: the same two 
dimers which are present after a plaquette flip on the corresponding kite has been 
performed, so as to increase the height function there. 

Any tiles meeting these criteria will work to reproduce perfect matchings on D n . The tiles 
we have chosen are the simplest possible ones: ali of the tiles are identical, ali of the faces 
are fiat, and the underside of each tile is planar. We have colored each tile according to its 
orientation. 

This construction appears to be highly artifìcial, so we should say a few words to explain 
why it is not. Our notion of height function comes from noncommutative Donaldson- 
Thomas theory [22|, in which we study the representation theory of a certain quotient of 
the path algebra CQ, the algebra of directed paths in Q (multiplication being given by 
concatenation when the paths share an endpoint, and the zero map when they do not). This 
theory was worked out first in the physics literature in [9 |. In particular, let Vccw and 
Vccw denote the set of counterclockwise (respectively, clockwise) directed paths around 
single faces in Q, and let W G Cq/[CQ, CQ] be the superpotential 

e n- e n- 

Let I w be the ideal generated by partial derivatives dW/da for each arrow a £ Q. The 
algebra that we study in noncommutative Donaldson-Thomas theory is CQ/Iw In this 
algebra, any two closed loops around a face with the same endpoint are identifìed; more- 
over, if p is a path from a to 6, and £ a and ij, are loops enclosing one face rooted at a and b 
respectively, one can show that 

(3) [4p] = [piai- 

(see jó'l for a proof). For our particular quiver Q, this algebra is non-commutative, but its 
center is described by a system of equations which cut out a cone on a del Pezzo surface. 
Our height change, as defìned, is designed so that any two paths which are equal in CQ/I W 
have the same height change. Indeed, consider the map 

*a : C Q /I W -> Q x Z> 

which sends a path to its endpoint in the universal cover Q together with its height change. 
We have chosen the defìnition of the height function h in such a way that ^> a is well- 
defìned, and in nice cases, namely when Q is consistent, is injective. Our quiver is in 
fact consistent. See fi] for a number of equivalent definitions and implications of consis- 
tency, as well as [3, 6 , 22] for an introduction to the portions of the representation theory 
of quivers relevant to this project. See also [7| for a physics perspective. 
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FIGURE 16. Weights on the lattice, before and after shuffiing. 
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8. The Generating Function 

We aim to give a generating function which counts, for each perfect matching, the 
number of tiles needed to construct M. Indeed, we can do slightly better: we can record 
some information about the orientations of these tiles. 

Theorem 7. 

n-l 

Z n (a,b,c) = Y[(l + a i b i c i+1 ) n - i (l + a i b i+1 c i+1 ) n - i . 



i=0 



n-l 

Ti— 3+1 



Z n+k (a, b,c) = n( 1+ Il ^ + a i+1 b i+1 éf-\ 

i=0 i=0 

Proof. Our strategy is to use a standard, folklore technique to encoding the weight of a 
perfect matching on D n as an edge weight function 

Wn-abc '■ E(D n ) — > {Laurent monomials in a, b, c} 

where a Laurent monomiai is a term of the form a % Wc k , i,j, k e Z. This was worked out 
explicitly for the square lattice in \24\. 

Consider a hexagon-shaped fundamental domain H of the dP 3 lattice, which has six 
squares, one square of each orientation. Assign weight 1 to both of the horizontal edges in 
H, as well as to the eight edges along the upper and lower boundary of H. Assign weight 
X to the easternmost vertical edge of the SE kite in H, and weight Y to the westernmost 
vellicai edge of the NW kite in H. Having done this, there is now a unique way to assign 
Laurent monomiai weights to the other edges in H such that a plaquette flip around one of 
the six faces / multiplies the weight of a perfect matching by the weight of /; it is shown 



in Figure 16 on the left. Denote these weig hts by a 



b,c- 

For n e |Z, let G be a minimal covering of D n by hexagons. Defìne w n[a b c {e) = 1 
for ali vertical edges e on the western boundary of G. This has the effect of setting Y = 
X(abc)^ 1 = 1 in iv^y bc ; for edges e in these hexagons, let 

W n ;a,b,c(e) = W abc,l;a,b,c- 

We can now defìne w n . a ^. c inductively for the remaining edges in G. Finally we restrict 
Wn-abc t° D n to obtain the desired weight function. 
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Then we have, for the NE-SW shuffle (which takes Z n to Z n+ i), 

Z n (a,b,c)=K ^ w n . a ^ c [Tx) eZ[a,b,c] 

7r perfect matching 



Now, when we perform domino shuffling on this weighted lattice, we allow each edge 



to retain its weight (see Figure 16 1. This can be achieved by a simple change of variables 
in the generating function: 

Z n+ i(a, (ac)- 1 ,(aò)- 1 ) = (1 + a) n+1 ■ KZ n (a, b, c) 

Here, K is a Laurent monomiai such that the Constant term of Z n (a, b,c) = 1. It is the 
inverse of the edge weight assigned to the minimum height function. 

Changing variables x = a,y — (ac) _1 ,z = (ab) -1 for the moment, and doing the 
similar computation for the NW — SE shuffle, we obtain the following straightforward 
system of recurrences: 



Z n+ i (x, y, z) = K'-{\ + x) n+1 ■ Z n (x, (zx)-\ (yx)' 1 ) 
Z n+l (x, y, z) = K" ■ (1 + z) n+1 Z n+ i ((yz)-\ (xz)-\z) 

for some Laurent monomials K', K". 

Let us solve these for the integer-order diamonds, leaving the half-integer order ones as 
a nearly identical exercise. We change variables again: x = a,y = b, z = c. It is also 
helpful to set q = abc, so that we obtain 

Z n+1 (a, b, c) = K'K"(1 + c) n+1 (1 + (òc)- 1 )^ 1 Z n ((bc)-\b, abc 2 ) 
= K'K"(1 + c)" +1 + ■ Z n (^Aqc) . 

thus, factoring out we have 




We know that the following are polynomials with Constant term 1 : 

Z n+l (a,b,c), (l + c)" +1 , + Z n (l,b,qc). 

Hence, the monomiai K'K" ( | j is equal to 1 and can be omitted from Q, giving 
the easy linear recurrence relation 

Z n+1 (a, b, c) = (1 + c) n+1 (l + Z n ^,6, qc\ , 

It is elementary to check that the first formula in the theorem statement satisfìes this first- 
order recurrence. As such, the theorem will follow once we check that Zi(a, b, c) is equal 
to (1 + c)(l + bc), which it is. 

□ 
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9. CONCLUSION 

We should reiterate that our main contribution is the domino shuffling algorithm itself, 
rather than our enumeration formulae. One other application of this algorithm is to a dif- 
ferent type of boundary conditions on this lattice, coming from Donaldson-Thomas theory, 
analogous to the pyramid parti tions studied in [24- 1 ; this example was mentioned in |[T7Ì . 
Indeed, we were originally inspired to study this problem because of its algebro-geometric 
connections. Accordingly, in a future paper, we will explain how to compute this generat- 
ing function. 

The reader may have noted that when n is an integer, the number of our order-n di- 
amonds is the square of the number of order n Aztec diamonds; one might hope for a 
bijective proof of this fact. However, this fact is not true of the generating functions, so 
such a proof will likely be very diffìcult to fìnd. 



Figure 17. A uniform random matching on an order 100 diamond 




We have used our shuffling algorithm to generate several perfect matchings of large 
diamonds, chosen uniformly at random from the set of ali matchings. An order- 100 is 



shown in Figure 17 One can clearly see the boundary approaching a smooth curve; indeed 
it is naturai 1 15 1 to guess that it is an algebraic curve, with boundary fluctuations given by 
the Airy kernel [12]. The simplicity of the domino shuffle that we have described suggests 
that the approach of [11] may be effective in calculating this limiting curve: the frozen 
boundary evolves according to a TASEP-like process which is likely amenable to analysis. 

The first author was supported in part by an undergraduate research scholarship from 
the Institut des Sciences Mathématiques (ISM). 
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